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THE LICK AGN MONITORING PROJECT: ALTERNATE ROUTES TO A BROAD-LINE REGION RADIUS 
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ABSTRACT 

It is now possible to estimate black hole masses across cosmic time, using broad emission lines in active 
galaxies. This technique informs our views of how galaxies and their central black holes coevolve. Unfortunately, 
there are many outstanding uncertainties associated with these "virial" mass estimates. One of these comes from 
using the accretion luminosity to infer a size for the broad-line region. Incorporating the new sample of low- 
luminosity active galaxies from our recent monitoring campaign at Lick Observatory, we recalibrate the radius- 
luminosity relation with tracers of the accretion luminosity other than the optical continuum. We find that the 
radius of the broad-line region scales as the square root of the X-ray and H/3 luminosities, in agreement with recent 
optical studies. On the other hand, the scaling appears to be marginally steeper with narrow-line luminosities. This 
is consistent with a previously observed decrease in the ratio of narrow-line to X-ray luminosity with increasing 
total luminosity. The radius of the broad-line region correlates most tightly with H/3 luminosity, while the X- 
ray and narrow-line relations both have comparable scatter of a factor of two. These correlations provide useful 
alternative virial BH masses in objects with no detectable optical/UV continuum emission, such as high-redshift 
galaxies with broad emission lines, radio-loud objects, or local active galaxies with galaxy-dominated continua. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1 . THE RADIUS-LUMINOSITY RELATION 

Over the past decade, interest in measuring supermassive 
black hole (BH) masses has intensified, as evidenc e mounts that 
BHs play a central role in galaxy evolution (e.g., ISilk & Reesl 
199% iHopkins et all l2006h . Locally, BH masses are mea- 
sured using stars, gas disks, or mega maser disks as dynam- 
ical tracers (e.g.. iGiiltekin et all 120091) . None of these tech- 
niques can currently reach beyond a few tens of Mpc. Thus, 
we resort to indirect mass estimates in actively accreting BHs 
to probe BH and galaxy coevolution at cosmological dis- 
tances. Studies of the BH mass and accretion-rate distri- 



butions both locally dGreene & Hoi 120071: ISchulze & Wis otzki 
20 1 0j) and at higher redshifts dWoo et alJl2006t iKollmeier et al. 
2006t lShen et alj2008rilWoo et al.l2008flVestergaard & Osmej 
2009; iKellv et alj 120091) . as weU as studies of possible 
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evolution in BH-bulge scaling relations (e.g., m 

20031: iTreu etal] | 2004t IWalter et al.ll2004t iPeng et alj|2 006a | 
Treu et al.ll2007h Isalviander et alj|2007t lAlexander et aljfe OQS 
Jahnke et al.ll2009b iGreene et al.ll2010l;lBennert et alj|2010l) , all 
rely on BH masses derived from active galactic nuclei (AGNs). 

BH masses derived from AGNs use the broad-line region 
(BLR) gas as the dynamical tracer, based on the assumption that 
the gas is primarily accelerated by the gravity of the BH. The 
gas velocity dispersion is derived from the broad-line width, but 
the BH mass estimate also requires the radius of the emitting re- 



gion. The best esti mate for its size comes from "reverberation" 
or echo mapping (Blandford & McKee 1982). Detailed spec- 
troscopic monitoring allows an estimate of the light-travel time 
through the BLR, by measuring the delay between variations in 
the continuum and line emission (see the recent compilation by 
Peterson et al. 20041) . This technique has a long history (e.g., 
Antonucci & Cohen] [1983T; IPeterson et al.l 119831: lUlrich et al.1 
19841 iGaskell & Sparkel 19861) , an d thus far has yielded reliable 
sizes for a few dozen sources (see lPeterson et al]|2004l) . 

Five reverberation-mapped sources show a 1 / y/R decline 
in velocity width ranging from C IV A 1549 to H/3, as ex- 
pecte d for a virialized BLR in a l/R potential (Kollatschnv 
120031: IPeterson et aT]|2004l) . Data from our Lick AGN Mon- 
itoring Project (LAMP), the subject of this paper, are consis- 
tent with the same assumption; when multiple Balmer lines 
are considered independently (e.g., Ha, H/3, Hy), all yield 
consistent estimat es for the so-called virial product, v 2 R/G 
dBentz et al.l 12010). On the other hand, other models, such 
as disk winds, would pr edict similar radial dependence (e.g., 
Murray & Chiang 1995). The importance of radiation pres- 
sure in supporting the BLR is currently a matter of de- 



bate as well (e.g iMarconi et al.1 120081 l2009t iNetzed [2009; 
iNetzer & Marzianij l2010t) . Despite these major uncertainties, 
the reverberation-derived BH masses correlate remarkably well 
with the lumi nosities and stellar velocity disper sions of their 
host bulges dBentz et al.ll2009bt IWoo et al.ll2010l) . In addition, 
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the very few existing direct dynamical measures of BH masses 
have so far turned out to be consistent with the reverberation- 
mapping virial estimates dDavies et alj2006tlOnken et alj2007t 
iHicks & Maikanll2008h . 

Since reverberation radii are usually not available, a sec- 
ondary estimate of BLR size is often obtained from the em- 
pirical correlation (the "radius-luminosity" relation) between 
AGN luminosity and BLR size, i? B LR oc iP dKaspi et alJ feOOO, 
20051 iBentz et al] ,2006. 2009a). With just a measurement of 
the AGN luminosity, typically L5100, and a broad-line width, 
typically FWHMh^, one can roughly estimate a BH mass as 
Mbh= fv 2 L,P/G. Here / is a scaling parameter that includes 
unknown information about the geometry and kinematics of the 
BLR. 

These so-called "single-epoch" virial BH masses are in- 
direct, and depend on a number of assumptions. Two em- 
pirically determined quantities fundamentally limit the accu- 
racy of the derived BH masses. One is /, which is cur- 
rently determined for ensembles of active galaxies through 
comparison between AGN-based masses and other estimates 
of M B h such as the M B h~ o* relation (e.g., [ Gebhardt et al 



uum luminosity. 



Table 1. Observations 



20001: iFerrarese et al.1 120011: iNelson et al.1 12004 lOnken et al 



2004; Gr eene & Holl200d: IShen et alJl2008at IWoo et alJl2010h . 
While there are good reasons to suspect that / may depend 
on physical prope rties of the BH such as accretion rate (e.g., 
ICollin et alJ .2006). reverberation-mapping campaigns have not 
yet succeeded in measuring / directly for individual objects. 
We are getting closer, however, since two-dimensional rever- 
beration mapping is growing more common and the velocity- 
resolved emission-line response strongly constrains / in in- 
dividual sources (e.g., iKollatschnvl [2003; Be ntz et alJ l2008t 
iDennev et al]|2009l) . 

The other empirically determined parameter is the slope of 
the radius-luminosity relation, /3, which is the subject of this 
paper. We are motivated to revisit this question thanks to our 
recent reverberation-mapping campaign, which has doubled the 
number of reverberation-mapped AGNs with Rblr ^10 light 
days. We do not consider the optical AGN continuum luminos- 
ity because the Hubble Space Telescope (HST) is required to 
spatially disentangle the AGN and galaxy continuum for these 
low-luminosity sources. The requisite HST imaging is under- 
way (GO- 11662, PI Bentz), and we will present the optical 
radius-luminosity relation in a future paper. Here we focus on 
other direct and indirect indicators of the AGN luminosity, in- 
cluding the X-ray luminosity and broad and narrow emission- 
line luminosities. 

There are practical reasons to consider other routes to de- 
termining BLR radii. For example, alternate relations are 
useful whenever the optical/UV continuum from the AGN is 
not measurable. This could occur when the AGNs are ra- 
dio loud, so that the optical/UV continuum is contaminated 
by synchrotron radiation, or when the galax y rather than the 
AGN dominates th e optical continuum (e.g., Wu et al. 2004; 
iGreene & Hol2005h . It has become common practice to use Ha 
or H(3 luminosities to calculate Rblr for high -redshift targets 
where the continuum is r arely detected (e.g.. lAlexander et al.l 
2008| IShapiro et al. 2009). Finally, remarkably, there is indi- 
rect evidence that broad-line widths measured from polarized 
line emissi on may provide a reasonable single- epoch virial BH 
mass (e.g., Zhan g^TaT] l200l iLiu et all 12009. C. Y. Kuo, in 
preparation). In these cases, hard X-rays or narrow emission 
lines are some of the only available proxies for AGN contin- 



Galax 


«BLR 


^•2-IOkeV 








Ref. 


. Dat _ 
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(3) 


(4) 
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fgj 


Mrkl42 


2.9± 1.3 


42.1 


41.64± 0.04 


41.2 


1 




5, 2007-11-24 


SBS1 1 16 


2.4± 1.3 




40.67± 0.04 


40.4 


1 






Arpl51 


4.1± 1.2 


42.5 


40.97± 0.07 


40.7 


1 




5, 2009-02-21 


Mrkl310 


3.7± 1.2 


42.3 


40.62± 0.04 


41.0 


1 




S, 2007-08-07 


Mrk202 


3.1± 1.6 


42.4 


40.48± 0.04 


40.4 


1 




A, 1999-11-09 


Mrk766 


6.2± 1.3 


42.9 


40.87± 0.02 


41.3 


1 


41.2 


S, 2006-12-28 


NGC4748 


5.6± 1.5 


41.3 


41 .01 ±0.02 


41.2 


1 




S, 2007-01-08 


NGC5548 


4.3± 1.3 


43.4 


41.35± 0.04 


41.6 


1 




5, 2007-6-19 


NGC6814 


6.7± 1.1 


42.2 


40.23± 0.04 


40.0 


1 




X, 2003-10-19 


Ivlrk335 


15.7± 1.3 


43.2 


42.14± 0.01 


41.5 


2 


41.0 




PG0026+129 


111.0± 1.3 


44.5 


42.92± 0.04 


42.9 


3 


42.2 




PG005 2+251 


89.8± 1.3 


44.8 


43.12± 0.05 


42.8 


3 






FairaU9 


17.4± 1.3 


44.0 


42.65 ± 0.02 


41.9 
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4 1 .5 




Mrk590 


25.6± 1.1 


43.6 


41.79± 0.05 


41.2 


2 


40.7 




3C120 


38.1± 1.6 


44.0 


42.36± 0.04 


41.9 


2 


42.5 




Arkl20 


39.7± 1.1 


43.8 


42.50± 0.03 


41.3 


2 






Mrk79 


15.2± 1.4 


43.5 


41.88± 0.03 


41.5 


2 






PG0804+761 


146.9± 1.1 


44.3 


43.26± 0.03 


42.3 


3 


41.7 




PG0844+349 


32.3 ± 1.6 


43.4 


42.52± 0.04 


41.6 


3 


41.2 




MrkllO 


25.5± 1.2 


43.9 


42.02± 0.02 


41.8 


2 






PG0953+414 


150. 1± 1.2 


44.6 


43.43± 0.02 


43.0 


3 


42.3 




NGC3227 


7.8± 3.8 


41.9 


40.49± 0.05 


40.4 


5 


40.3 




NGC3516 


6.7±2.2 


43.1 


41.33± 0.05 


40.8 


5 


41.0 




NGC3783 


10.2± 1.3 


43.0 


41.49± 0.03 


41.2 


6 


40.9 




NGC405 1 


5.8± 1.4 


41.4 


39.70± 0.08 


39.7 


7 


39.6 




NGC4151 


6.6± 1.2 


42.7 


4U0± 0.09 


41.6 


8 


40.8 




PG121 1+143 


93.8± 1.5 


43.7 


42.99± 0.05 


42.4 


3 


41.4 




PG 1226+023 


306.8± 1.3 


45.8 


44.10± 0.03 


43.0 


3 


42.8 




PG 1229+204 


37.8± 1.7 


43.4 


42.31 ±0.05 


41.9 


3 


41.4 




NGC4593 


3.7± 1.2 


42.8 


40.89± 0.11 


40.6 


9 


40.4 




PG 1307+085 


105.6± 1.5 


44.4 


43.14± 0.05 


42.8 


3 


42.1 




Mrk279 


16.7± 1.3 


43.7 


41.96± 0.04 


41.5 


10 


41.4 




PG141 1+442 


124.3± 1.7 


43.2 


42.83± 0.02 


42.0 


3 


41.6 




NGC5548 


18.0± 1.0 


43.5 


41.73±0.05 


41.6 


11 


40.9 




PG 1426+0 15 


95.0± 1.5 


44.0 


42.80± 0.04 


42.2 


3 


41.9 




Mrk817 


21.8± 1.1 




41.97± 0.03 




2 






PG1613+658 


40. 1± 1.5 


44.4 


42.9S± 0.03 


42.4 


3 


42.4 




PG1617+175 


71.5± 1.6 




42.7 1± 0.05 




3 






PG 1700+5 18 


251. 8± 1.2 




43.74± 0.02 




3 






3C390.3 


23.6± 1.3 


44.2 


42.29± 0.03 




12 






Mrk509 


79.6± 1.1 


44.1 


42.60± 0.00 


42.3 


2 


41.9 




PG2 130+099 


22.9± 1.2 


43.7 


42.74± 0.04 


42.1 


3 


42.0 




NGC7469 


4.5± 1.2 


43.2 


41.76± 0.02 


41.7 


13 


41.4 





Note. — Col. (1): Galaxy name. Note lhai references lor all data sels arc included in §3. Col. (2): BLR 
radios (liglrl days). Col. (3): X-ray luminosity (erg s~'). Errors arc not included, since they tire dominated 
by variability, but we assume 2(l'/i uncertainties on all values (sec §3.1). Col. (4): H.fcl luminosity (erg s _1 ). 
Col. (5): [O lit] A5007 luminosity (erg s H ); assumed errors arc 0.15 dcx, based on differences in literature 
measurements (§3). Col. (6): References for the [O tit] fluxes: (1) Bentz et al. (2009); (2) Peterson et al. 
(1998); (3) Boroson & Green (1992); (4) Winge et al. (1996); (5) Denney et al. (2010); (6) Stirpe et al. 
(1994); (7) Peterson el al. (2000); (8) Kaspi et al. (1996); (9) Dietrich el al. (1994); (10) Santos-Lied et al. 
(2001); (11) Peterson el al. (1991); (12) Dietrich et al. (1998); (13) Collier et al. (1998). Col. (7): [O 
tv] A25.8 /rm luminosity (erg s _1 ); assumed errors arc 0.15 dex as for [O ill] above. Col. (8): Instrument for 
X-ray luminosity: .V -Swift. C -Chimdm. X =XMM-Nm , nm. A =ASCA. Below the line all X-ray luminosities 
are taken from Kaspi ct al. (2005). l or the LAMP objects we also include the dale of observation. 



1.1. Which Luminosity Best Predicts BLR Size ? 

If the spectral energy distributions (SEDs) and the den- 
sity distribution in the BLR are independent of luminosity, 
then we expect that the BLR size will scale simply with the 
square-root o f the photoionizing luminosity, 7?blr oc \/L (e.g., 
Netzer 1990). The most recent calibrations of the radius- 
luminosit y relation have all been consistent with this simple 
relation dBentz et al.ll2006ll2009al) . 

If the SEDs were really independent of luminosity, the pho- 
toionizing luminosity could be estimated from a measure of 
the AGN continuum at almost any wavelength. However, 
there are observational indications of luminosity dependence 
in SEDs. Indeed, on theoretical grounds we might also ex- 
pect smaller and hotter accretion disks arou nd lower-mass BHs 
(e.g.. lShieldslll978HZheng & Malkanlll993l) . While the equiva- 
lent widt h of H/3 is constant in h igh-lumi nosity active galax- 
ies (e.g.. ISearle & Sargentl 119681) . both ICroom et all d2002l) 
and IGreene & Hoi ( 20051) see evidence for a weak inverse- 
Baldwin effect in H/3 at low luminosity. Furthermore, the in- 
crease in q ox with UV luminosity (e.g.. [Avni & Ta nanbaum 
119821: ISteffen et al.l 120061: iDesroches et al.1 120091) suggests 
luminosity-dependent changes in the SED. Finally, the rel- 
ative strengths of the "big blue bump" and the X-rays de- 
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pend on Eddington ratio, with the latter dominating at lower 
Lhn. lLv M dMalkan & Sargent! 19821: 1 Vasudevan & Fabianl2007t 
IVasudevan et alj|2009j). For a review of the situation at yet 
lower Lboi/^Edd see iHol 120081 Therefore, in this paper we 
will consider several observables which may correlate with, 
and thus be used to estimate, the ionizing luminosity. Specifi- 
cally, we consider the following proxies for the AGN luminos- 
ity: hard X-ray luminosity (Lz-iq kev)> H/J luminosity (Z-h/3), 
narrow [O III] A5007 A luminosity (L[o nrj), and narrow [O IV] 
A25.8 /im luminosity (L[o iv])- 

Throughout we assume the following cosmological param- 
eters to calculate distances: Ho = 100/z = 70 km s" 1 Mpc -1 , 
fl m = 0.30, and n A = 0.70. 

2. THE LICK AGN MONITORING PROJECT 

The new measurements that motivate this work result from 
LAMP, a dedicated monitoring campaign of 13 AGNs (includ- 
ing the well-studied Seyfert galaxy NGC 5548). We specif- 
ically focused on nearby (redshift z < 0.05) Seyfert galaxies 
with low luminosities (AL5100 < 10 43 erg s" 1 ) and probable 
BH masses in the range 10 6 — 3 x 10 7 M , since this luminos- 
ity and mass regime had not been explored fully in the past. 
Spectroscopic monitoring was carried out with the Lick Ob- 
servatory 3 m Shane telescop e over a nearly contiguous 64 
day period dBentz et alj|2009a) . while photometric monitoring 
was per formed over a long er period utilizing four smaller tele- 
scopes (Wals h et al.l 12009ft . We successfully measu red BLR 
radii based on H/3 for nine objects dBentz et alj2009cj). reporte d 
lags in multiple other Balmer transitions (Bentzetal. 2010), 
and succeede d in measuring velocity-r esolved lags in at least 
three sources dBentz et alj|2008l l2009d) . Finally, we revisited 
the calibration of reverberation-m apped BH masses using the 
Mbu-<t* relation dWoo et al.ll2.Q10 ). For the purpose of this pa- 
per, we focus on BLR radii based exclusively on H/3 lag times. 

3. LUMINOSITIES AND BLR RADII 

The BLR light-c rossi ng times used here are presented by 
iBentz et all d2009al) and iBentz et all d2009d) for the previous 
reverberation-mapped and LAMP AGNs, respectively. We note 
that improved lag measur ements were more re cently reported 
for a subset of galaxies bv lDennev et alj d2010h . We have con- 
firmed that the radius-luminosity relation based on H/3 does not 
change with the inclusion of their lag values, but continue to use 
the old measurements f or temporal consiste ncy with the X-ray 
obser vations. We follow lBentz et alJ d2009a) and lPeterson et al.l 
(2004) and remove IC 4329A from the sample due to uncertain- 
ties in the measurements. Throughout we will refer to the sam- 
ple of active galaxies with reverberation mapping, excluding 
the LAMP targets, as the "non-LAMP" objects. We describe 
the origin of the AGN luminosities in this section (Table 1). It 
is useful to remember that "BLR size" here actually refers to 
the time of peak response of the H/3 emitting gas relative to the 
continuum, multiplied by the speed of light. Had the experi- 
ment been done with C IV A 1549, for instance, the sizes would 
have been smaller, but the widths larger. When calculating ef- 
fective BH masses, it is important to match the species used to 
measure velocity dispersion with the radius relation calibrated 
for the same species. 

Of all the luminosities we discuss, only the broad 
H/3 line luminosity (Lh/3) is measured as part of the 
reverberation-mapping campaign, simultaneously with 
the radius measurement. It should provide a fairly 



direct and unbias ed probe o f the p hotoio nizing con- 
tinuum . Both IBentz et alj d2009d) and iKaspi et all 
J2001 tabulate average L H/ 9 measured in the same way. 
Table 2. Fits 



L 


a 


P 


eo 


Fit 


(1) 


(2) 


(3) 


(4) 


(5) 


i2-10kev/10 43 


0.09±0.05 


0.52±0.05 


0.26 


C 


^2-lOkev/lO 43 


0.09±0.06 


0.52±0.05 


0.25±0.05 


ML 


Wio 43 


O.85±0.O5 


0.53±0.04 


0.22 


C 


iH/3/10 43 


0.85±0.06 


0.53±0.04 


0.22±0.04 


ML 




0.86±0.06 


0.53±0.04 


0.22 


C 


Wio 43 


0.86±0.06 


0.53±0.04 


0.22±0.04 


ML 


i[oni]/10 42 


0.53±0.06 


0.62±0.07 


0.29 


C 


i[om]/10 42 


0.52±0.06 


0.61±0.07 


0.30±0.05 


ML 


^[oiv]/10 42 


0.76±0.10 


0.58±0.11 


0.35 


C 


^[Oiv]/10 42 


0.75±0.10 


0.58±0.11 


0.35±0.07 


ML 



Note. — Fits to log (^blr/ lOpc) = a + /3 logL, for each luminos- 
ity. 

Col. (1): Luminosity measure (erg s _1 ). L H ^ b /10 43 has the narrow 
H/3 emission removed. Col. (2): a. Col. (3): 0, Col. (4): Intrin- 
sic scatter. Maximum-likelihood fits are those with error bars on the 
intrinsic scatter. Col. (5): Fit type, either C=\ 2 fit or ML=maximum- 
likelihood. 

Note that we present results based both on the total H/3 lumi- 
nosity (narrow and broad combined) and the broad H/3 luminos- 
ity alone. The results are basically identical, since the median 
luminosity difference is less than 5%. Although it would be 
useful to examine Ha as well, uniform measurements do not 
exist for the non-LAMP sample, and thus we must await future 
work. 

The [O III] luminosities for the LAMP sample itse lf are mea- 
sured fr om the Shane spectra and are presented bv lBentz et al.l 
d2009cT) . For the non-LAMP targets, we draw from previ- 
ous reverberation mapping campaigns for the local galaxies. 
For the more distant an d luminous Palomar-Green quasars 
dSchmidt & Greenl 1 1 983b . we combine the equivalent- width 
measurements of [O III] from iBoroson & GreerJ (Q992) with 
the continuum flu xes o f iKell ermann et al. (1989) as given by 
iHo& Pengl d2001l) and lGreene et al.l d2006l) . Table 1 contains 
all measurements, including relevant references. Note that the 
[O III] luminosities have not been corrected for extinction. The 
formal errors for the [O III] measurements are in the range 
2-15%. However, we find a median difference of ~ 40% be- 
tween differen t literature valu es (relying predominantly on the 
compilation of I Whitflel [r992) . The values used here, from pre- 
vious reverberation mapping campaigns, are typically smaller 
than those compiled by Whittle. Thus, in our fitting we adopt an 
uncertainty of 0.15 dex in the line luminosities, as an estimate 
of the impact of various systematic effects discussed below. 

A large fraction of the non-LAMP sources have [O IV] A25.8 
^m luminosities available in the literature from the Spitzer 
Space Tel escope. We draw p redominantly from the measure- 
ments of Das vra et alj d2008l) . which cover a large fraction of 
the non-LAMP reverberation-mapped sample. We take the 
measurements for Mrk 766 and Mrk 335 from the work of 
Tommasin et al. (2010) and that of NGC 3516 from the work 
of iGallimore et al. ( 20101) . With the exception of this last, 
all were taken with the high-resolution grating. A similar 
comparative exercise as above, this time with the compila- 
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tions of Tommasin et al. | d2010t) . iGallimoreetafl d2010l) . and 
IVeilleux et al.l ( 20091) . yields a lower uncertainty estimate (~ 
10%) than for the [O III] lines. Presumably greater agreement 
is reached because the data sets were in many cases identical, 
and so we use a value of 40% as above. 

The X-ray luminosities for the non-LAMP sourc es are taken 
directly from the compilation of Kasp i et ail d2005l) . They are 
derived from a variety of literature sources, but generally are 
based on Ad vanced Satellite for Cosmology and Astrophysics 
(ASCA) data ( iTanakaet al.ll 1994 . While many of these targets 
have more recent XMM-Newton or Chandra observations avail- 
able in the archive, the ASCA measurements are actually closer 
in time to the reverberation-mapping campaign. Thus, we adopt 
the X-ray luminosities presented by Kaspi et al. in all cases. 

The X-ray luminosities for the LAMP objects come from a 
range of sources. In all cases we adopt the observation clos- 
est in time to that of our campaign (Spring 2008). Six objects 
(Mrk 142, Apr 151, Mrk 766, Mrk 1310, NGC 5548, and NGC 
4748) have X-ray luminosities from Swift. All were observed 
between 2007 and 2009. The Swift data were extracted using 
the xselect task as part of the HEASARC tool-set. Each source 
was extracted from a circular region with a radius of 20 pixels 
(~ 47"). Background rates were negligible in all cases. Count 
rates were converted to fluxes assuming a power-law spectrum 
with r = 1.8 (E oc E~ r ) and no internal absorption. In the few 
cases with multiple epochs (e.g., Mrk 766 and NGC 5548), we 
analyze the longest observation where the galaxy center is close 
to the image center. 

The remainder of the LAMP AGNs only have heterogeneous 
measurements available in the literature. In on e case (NGC 
6814) we use the XMM-Newton slew survey (ISaxton et al.l 
l2Q08h and in anoth er (Mrk 202) we resort to an ASCA obser- 
vation from 1999 dUeda et al]|2005l) . where aperture photome- 
try yields a count rate that is converted to a flux assuming only 
Galactic extinction and V = 1 .7. The luminosity for NGC 68 14, 
from the slew survey, was derived in a similar fashion, using 
the same spectral model. The only difference is that the flux is 
reported for 2-12 keV. We use we bPIMMSH to calculate the 
2-10 keV flux assuming our spectral model. 

With many years of comprehensive monitoring, NGC 5548 
is a special case and warr ants extra atten t ion. T here are fourteen 
non-L AMP epochs from lPeterson et all d2002l) and lBentz et al.l 
d2007l) . In Figure lb we show all fifteen epochs in gray for ref- 
erence. Currently, NGC 5548 is in a very low luminosity state. 
The LAMP measurement differs by a fact or of ~ 4 from th e 
weighted average of all other epochs (e.g., Bentz et al. 2009c). 
Unfortunately, we have only two epochs of X-ray data for this 
source, and, given narrow-line region sizes of hundreds of pc, 
the [O III] luminosity is presumed constant over timescales of 
months. For the purposes of fitting, we adopt the wei ghted av- 
erage lag, 18 ± 0.6 light-days, from lBentz et al.l (l2009al) as the 
non-LAMP point. The early X-ray data are from ASCA and 
were taken in 1993, when the lag was measured to be 13+j \ 
light-days dPeterson et al.ll2002l) . If we rather adopt the latter 
value in our fitting of the X-ray radius-luminosity relation (§4), 
it makes no difference to our results. Since the BLR size of 
NGC 5548 has been observed to change on timescales short 
compared to changes in the narrow-line flux, it is also interest- 
ing to note that there is a scatter of 0.2 ± 0. 1 dex in the logarithm 

'http://heasarc.gsfc.nasa.gov/ 
10 http://heas arc. gsfc.nasa.gov/Tools/w3pimms. html. 
1 1 http://tartarus.gsfc.nasa.gov/. 



of the ratio of lag to [O III] luminosity across the fifteen epochs. 

3.1. Systematics: X-ray Variability, Aperture Effects, and 
Extinction 

Each of the luminosities we consider comes with its own 
complications. In the case of the narrow emission lines, they 
have been photoionized by the average continuum luminosity 
over the past ^100 years, during which time Rblr may vary 
significantly. On the other hand, the X-ray emission region is 
more compact than the optical emitting region, and thus varies 
on shorter timescales than changes in Rslr occur. With nonsi- 
multaneous observations, we may introduce significant scatter 
into the /^blr-^-io keV relation. In addition, X-ray variability 
timescales depend systematically on Mbh and luminosity (e.g ., 
lO'Neill et al]|2005b iMcHardv et aTll2006t iMiniutti et ail l2009). 
Thus, it is at least conceivable that some systematic bias is in- 
troduced into the /?blr-^2-io keV relation. We investigate that 
possibility here. 

We start by considering all multi-epoch data available for 
reverberation-mapped sources from the Tartarus databasJH. 
The benefit of Tartarus is that the fluxes have been derived 
from the ASCA X-ray spectra in a uniform way. Spectral fits 
to the hard X-rays (2-10 keV) are performed, with the region 
around Fe Ka masked and including possible internal absorp- 
tion (which is small in this spectral region). Ten of the non- 
LAMP targets have multiple epochs of observations spanning 
more than one year in the Tartarus database. They include 3C 
120, Fairall 9, Mrk 509, NGC 3227, NGC 3516, NGC 3783, 
NGC 4051, NGC 4151, NGC 4593, NGC 5548, NGC 4269, 
and PG 1226+023. The typical cadence is a few observations 
per year. 

For each object we calculate a mean and root-mean square 
(rms) flux using the Tartarus database. We find variability am- 
plitudes of 5-80% (one standard deviation) over the 1-7 yr 
timescales probed by these observations. NGC 3516 is the 
target with the highest variability amplitude (80%). The ma- 
jority of objects do not vary even by a factor of two on these 
timescales. The median amplitude of variability is ~ 20%. 
Thus, the level of intrinsic variability in the X-ray luminosity of 
most Seyferts is usually too small to impact the Rblr-Li-w keV 
relation. For fitting purposes we thus adopt 20% uncertainties 
in all X-ray fluxes. 

We perform a second check using artificial light curves. Our 
goal is to investigate whether systematic changes in break time- 
scale will lead to a bias in o ur derived ^BLR-^2-iokeV relation. 
We use the prescriptions of iTimmer & Koenigl d 1995b to gen- 
erate mock light curves with an input power spectrum of vari- 
ability. We generate a family of light curves, each of 5 yr dura- 
tion, and each with a characteristic break in the power spectral 
density function. The break timescales range from 0.01 to 30 
days, which is similar to the range of 0.005 to 30 days seen i n 
reverberation-mapped sources (e.g., lUttlev & McHardv1l2005h . 
For simplicity, all light curves have a power-law slope of a = -2 
(P oc f a ) at frequencies above the break frequency and a slope 
of (3 = -1 (P oc f°) at frequencies below the break. Each arti- 
ficial light curve is "observed" 2000 times with 10 ks duration. 
The signal-to-noise ratio is taken to be 100 (but does not im- 
pact the results) and the assumed variability amplitude is 20% 
to match the Tartarus average above. We then look at the spread 




40 41 42 43 40 41 42 

lo S ( L [o in] / er g s_1 ) lo g ( L [o iv] / er g s_1 ) 

FIG. 1. — (a) Fit to the L2-10 keV _ ^BLR relation, including both LAMP sources (open circles) and non-LAMP sources (filled dots). Our maximum-likelihood 
fit (solid; = 0.52 ± 0.05) and the fit of Kaspi et al. (2005; dashed) agree in this case. The intrinsic scatter (Table 2) is ~ 0.3 dex. NGC 5548 is included twice, 
and indicated with a red box in all panels, (b) As in (a), but here using the H/3 rather than X-ray luminosity. I n this case, our max imum-likelihood fit (solid; 
/3 = 0.53 ± 0.04) is significantly shallower than that of Kaspi et al., bringing the slope into agreement with that of Bentz et al. 1 2009a) for the optical continuum. 
For reference, we show all fifteen epochs of monitoring for NGC 5548 in gray, (c) The Kblr-£[0 III] relation. As above, LAMP sources are open circles while the 
non-LAMP sources are small black circles. The solid line is our best maximum-likelihood fit (/3 = 0.60 ±0.07). (d) As in (c), but here using the [O iv] 25.8 fim lumi- 
nosity rather than the [O 111] luminosity. The arrow indicates the upper limit on the [O iv] luminosity of PG2130+099. Our maximum-likelihood fit (J3 = 0.69±0.13) 
is shown as a solid line. 
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in derived X-ray fluxes as a function of break timescale. There 
is no change in the width of the distribution of mean fluxes for 
breaks ranging from 0.01 to 30 days. Therefore, we do not 
expect any systematic errors in our X-ray flux estimates as a 
function of mass or luminosity based on trends between bre ak 
timescale and mass or Eddington ratio dMcHardv et alj|2006l) . 

It is also worth discussing the primary sources of systematic 
errors in the fluxes of narrow emission lines. First, aperture cor- 
rections can be significant, since narrow-line regions (NLRs) 
have sizes o f hundreds of parsecs and are often spatially re- 
solved ( e.g.. IWhittldll992t iBennert etafl l2002t ISchmitt et all 
2003a b; Greene et al. 2009). Objects that are closer are more 
susceptible to aperture effects. The LAMP targets were ob- 
served with a slit of width 4", while the non-LAMP targets 
generally come from apertur es of width 2-20" (altho ugh most 
are larger than 4"), and the lBoroson & Greenl (Tl992) observa- 
tions were taken wit h a 1 "5 slit. Between the c ompilations of 
IBennert et all d2002l) and ISchmitt et ail d2003al) . we find NLR 
sizes for four of the PG quasars, two LAMP objects (NGC 5548 
and Mrk 766), and five other non-LAMP objects (Mrk 590, 
NGC 35 16, Mrk 79, NGC 3783, NGC 4593). The NLR sizes of 
these galaxies are 0.1-2", with a median size of 0"3. The PG 
quasars all have NLR sizes < 1 "2. Nominally, based on the ob- 
jects with measured NLR sizes, we expect minimal loss of light 
due to aperture effects. Observing c onditions will also l ead to 
some slit losses, in the case of the iBoroson & Greenl (119921) 
observations, but not at a level that is significant compared to 
internal extinction (see below). To address the rest of the sam- 
ple, we estimate NL R sizes using the size-luminosity relation of 
ISchmitt et ail d2003bl) . We find that the typical expected NLR 
size (for both LAMP and non-LAMP sources) is — 1"3. The 
PG quasars, which were observed with the smallest slit, all have 
expected sizes in the range 0"5- 1". We are in even less danger 
with the L[o ivj measurements, given th eir larger aperture size s 
(5-20") and compact emission regions dMelendez et al.ll2008l) . 

The next important concern is internal extinction. It has 
been shown many times, particularly for obscured AGNs, that 
the [ O III] luminosity can be significantly extinguished by dust 
(e.g.. lMalkanlll983l; iMulchaev et all [19941 iNetzer et al.i r2006; 
Melende z et alj|2008l) . On the other hand, L[o ivj is relatively 
insensitive to extin ction, making it a h i gher fidelity luminos- 
ity in dicator (e.g jMelendez et al]|2008l : iDiamond-Stanic et alj 
120091) : if the L[o mj measurements are compromised, we still 
expect to find reasonable results for Ljoiv]- Finally, there is 
the possibility of contamination from star formation. In prin- 
ciple this is possible for L[o nrj, but note that the [O III]/H/3 
intensity ratio is considerably lower in high- metallicity star- 
formi ng galaxies than in active galaxie s (e.g., iBaldwin et al.l 
119811) . As shown by iKauffmann et alj d2003l) . the expected 
level of L[o nrj contamination from star formation in local ob- 
scured active galaxies is low (< 10%). 



4. FITS 

Our primary goal is to calibrate the relation between BLR 
size and various indicators of nonstellar (AGN) luminosity. We 
fit to the standard relation (/?BLR/101t-days) = a+/?logL, where 
L here is derived from H/3, L2-iokeV, and narrow emission-line 
luminosities. For comparison with recent literature, we utilize 
two primary fitting schemes. The fir st is a y 2 minimizatio n 
technique similar to that presented by Tremai ne et alj d2002l) . 



The following \ 2 function is minimized: 



XT 



E 



(Ri-a-pLj) 2 
4, + /3 2 4 



(1) 



Intrinsic scatter is accounted for by replacing e«, with e = 
(e^ + e 2 ,) 1 ' 2 , where eo (the intrinsic scatter) is chosen such that 

X?=l- 

In addition, we use a m aximum-likelihood technique adapted 
from Gul tekin et al.l d2009l) . For simplicity we assume that both 
the measurement errors and the intrinsic scatter have Gaussian 
distributions. For a set of observed points (/?,,L,), we maximize 
the total likelihood, 



C = Y[li(Ri,Li). 



(2) 



In the presence of measurement errors, if the likelihood of mea- 
suring a BLR radius R, for a true radius R is Qi(Rj\R)dRj, and 
the probability to have a true radius R given L, is P, then for a 
given observation the likelihood is 



k = J Qi{Ri\R)P{R\LddL. 



(3) 



We assume that both Q and P have a log-normal form. Un- 
certainties in the independent variable (luminosity) are derived 
from Monte Carlo simulations and are always small. Fits using 
both methods are given in Table 2, with the first line showing 
the x 2 method. In all cases the results of the two fitting methods 
are indistinguishable. 

It is interesting to note that the X-ray and H/3 relations 
are now consisten t with a slope of Rblr oc \fL. In contrast, 
Kas pi et al.l ((2005) report a slope of 0.7 for the X-ray relation 



(R 



BLR 



0.7 

2-10keV 



). We should note, however, that when they fit 



an average lag for each object and used only H/3 lags (the most 
directly comparable case to what we have done here), they find 
a slope of 0.5 (Rblr oc Z^-iokev)- Their reported Rblr- Lap 
slope is steeper, Rblr oc L^ 9 . With our improved data, we find 
that both relations are consistent with a slope of 0.5. Thus, 
the simplest assumption, that AGN SEDs and BLR densities 
are independent of luminosity, appears to apply, at least for 
the present sample and to the level of precision that can be 
tested by our data. One goal of ongoing reverberation-mapping 
campaigns should be to investigate whether there are physical 
regimes (e.g., in BH ma ss or luminosity) for which this assump- 
tion does not hold (e.g jGreene & Holl2009l) . 

There is tantalizing evidence, in contrast, that the narrow-line 
relations may have a steeper slope, although with low signifi- 
cance. Here we explore possible interpretations of this result, 
should it turn out to be significant. Above we discussed vari- 
ous sources of contamination of the NLR luminosity, including 
redshift-dependent aperture correction, extinction, and star for- 
mation. Aperture effects go in the wrong direction to explain 
the steeper slope, while extinction seems implausible because it 
would have to impact the L[o rvj measurements as strongly than 
the L[o mj measurements, contrary to normal reddening laws. 
Star formation could artificially boost the NLR luminosities at 
the low end. However, we do not believe the L[o nrj contami- 
nation could be more than ~ 10% on average, while the values 
need to be boosted by factors of 2-3 to impact the slope on a 
logarithmic scale. Therefore, the steeper slope, if real, is more 
likely explained by physical effects rather than measurement 
errors. It would most naturally arise from the measured lumi- 
nosity dependence in the relation between NLR and bolometric 
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luminosity. We are not the first to report this trend. For in- 
stance, iNetzer et al.l d2006l) find that L [0 m oc I%!°^v> while 
IMelendez et all d2008l) find L [0 IV ] oc ^-lOke'v and ^[o ni] °c 
^[OW] 1 - ^ ^blroc ^2-10 kev°' 5 > th en based on Netzer et al. we 
would expect /?blroc £[0 ni]°' 7±<U , which is consistent with our 
finding. The slope we measure i n the Rblr-Liq ivi rel ation is 
also consistent with the results of lMelendez et al.l d2008l) . Thus, 
the possibility of a steeper slope is plausible. For some reason, 
quasars are less efficient at powering an NLR than are the less 
luminous Seyfert nuclei. 

There is now compelling evidence that bolometric correc- 
tions depend on the Eddington ratio Lboi/^Edd, where the Ed- 
dingt on luminosity for 1 is taken to be 1.25 x 10 38 e rg s -1 
(e.g., IVasudevan & Fabianll2007l: IVasudevan et alj|2009t) . It is 
worth seeking correlations between radius-luminosity relation 
residuals and the Eddington ratio. BH mass measurements are 
provided by the reverberation-mapping campaigns. In the case 
of the non-LA MP sources, w e take the Eddington ratios from 
the study of Vasu devan et alj d2009l) . who have measured si- 
multaneous SEDs ranging from the optical to the X-ray using 
XMM-Newton. We do not yet have full SEDs for the LAMP 
sample, and so we use a single-band observation and a bolo- 
metric corre ction. We adopt £2-10 ke v and the bolometric cor- 
rection from IVasudevan et ail d2009ft . The bolometric correc- 
tion depends on the Eddington ratio and we assume a value of 
30, as appropriate for sources with Lboi/£Edd~ 10%. The result- 
ing Eddington ratios are in the range 0.001-1, but are strongly 
peaked at ~ 0.1. 

We seek correlations between the residuals in Rblr around 
the mean Rblr-L relations and the Eddington ratio. The non- 
parametric Kendall's r is calculated (within IRAtf^l) for all re- 
lations. In no case do we find evidence for a correlation be- 
tween the Rblr-L residuals and Lboi/^Edd- The probability of 
no correlation is in the range P = 0.3-0.8. 

Although we do not know its origin, it is interesting to ex- 
amine the intrinsic scatter for each fit. As expected, the intrin- 
sic scatter is lowest when Lup is used, presumably because of 
both temporal and spatial coincidence. On the other hand, the 
relations based on both the X-rays and narrow emission lines 
have comparable scatter. One might expect higher scatter in the 
narrow emission-line relation due to the unquantified role of in- 
ternal extinction and aperture effects. Furthermore, the narrow- 
line emission cannot respond to changes in accretion luminosity 
on timescales of a month, while we know that Rblr does. Note 
that NGC 5548 has shown Rblr variability at the factor of four 
level, and yet the overall relations only have an intrinsic scatter 
of a factor of two. Once the L5100 measurements are in hand, 
it will be interesting to see whether the intrinsic scatter is min- 
imized using the optical continuum luminosity or, indeed, the 
bolometric luminosity. 

5. SUMMARY 

We explore radius-luminosity relations based on AGN lumi- 
nosities other than the optical continuum. The time is right 
to revisit these relations because of a new sample of low- 
luminosity, low-mas s AGNs with rever beration mapping from 
the LAMP project dBentz et alj|2009d) . We consider X-ray, 



broad H/3, narrow [O III], and narrow [O IV] luminosities. 
These relations are designed for use in estimating BH masses 
when optical continuum luminosities are not available. Rel- 
evant situations include local AGNs with galaxy-dominated 
spectra, and possibly radio-loud objects, various high-redshift 
active galaxy populations (such as submillimeter galaxies), and 
heavily obscured AGNs with detected broad polarized emis- 
sion. Furthermore, any differences in slope or intrinsic scatter 
between relations based on different luminosities may indicate 
SED differences in AGNs as a function of luminosity or BH 
mass. 

We find that the Rblr-Li-io kev and Rblr-Lhp relations are 
well fit with a slope of /?blroc \fh. This is the slope expected 
if AGN SEDs and BLR densities are independent of luminos- 
ity. On the other hand, the narrow emission lines show tenta- 
tive evidence for a steeper relation, /?blroc L 6 . Intriguingly, 
these slopes are consistent with previous results showing that 
L[o r m/Lx and Z40 ivi/^x decreases with incre asing luminosity 
(e.g.. [Netzer et aljl2006t IMelendez et al.ll2008l) . We find no ev- 
idence for a correlation between Rblr-L residuals and Edding- 
ton ratio. In fact, the intrinsic scatter in all relations is surpris- 
ingly small. On the one hand, the X-rays are variable on short 
timescales, but, as we show, that does not translate into signif- 
icant errors in the /?blr-^2-io kev relation. On the other hand, 
the narrow emission-line luminosities do not respond at all to 
state changes on timescales of a year. Thus, we find it surpris- 
ing that even in these cases the intrinsic scatter is only at the 
factor of two level. Still, this scatter translates directly into 
uncertainties in the BH masses (e.g., IVestergaard & Peterson! 
2006J iMcGill et al.ll2008h . As the reverberation-mapped sam- 
ples increase, it should become possible to search for evidence 
of secondary parameters that might allow one to decrease the 
total scatter, thereby increasing the fidelity of our BH mass es- 
timates. 
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